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Abstract: We consider the renormalization group improvement in the theory of the Stan- 
dard Model Higgs boson playing the role of an inflaton with a strong non-minimal coupling 
to gravity. It suggests the range of the Higgs mass 135.6 GeV < Mh ^ 184.5 GeV com- 
patible with the current CMB data (the lower WMAP bound on n^), which is close to the 
widely accepted range dictated by the electroweak vacuum stability and perturbation the- 
ory bounds. We find the phenomenon of asymptotic freedom induced by this non- minimal 
curvature coupling, which brings the theory to the weak coupling domain everywhere ex- 
cept at the lower and upper boundary of this range. The renormalization group running 
of the basic quantity Aj - the anomalous scaling in the non-minimally coupled Standard 
Model, which analytically determines all characteristics of the CMB spectrum - brings A/ 
to small negative values at the inflation scale. This property is crucial for the above results 
and may also underlie the formation of initial conditions for the inflationary dynamics in 
quantum cosmology. 
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1. Introduction 

The task of constructing a fundamental particle model accounting for an inflationary sce- 
nario and its observationally consistent predictions for the CMB was undertaken in Q 
and, more recently, in the series of papers g |, |,|, |. While the relatively old work 
||T| had suggested that due to quantum effects inflation depends not only on the inflaton- 
graviton sector of the system but rather is strongly effected by its GUT contents, the paper 
by Bezrukov and Shaposhnikov [Q] transcended this idea from the remote land of Grand 
Unification Theory to the seemingly firm ground of the Standard Model (SM) with the 
Higgs field playing the role of an inflaton. This has provoked new interest in a once rather 
popular [0, ^, |lO|, ID, |ll|, ^] but then nearly forgotten model with the Lagrangian of the 
graviton-inflaton sector given by 

LiOf.u,'!^) = ^ (Ml + m') R - ^IV^P - Vm), (1-1) 

y(|$|) = ^(|^|2 _^2^2^ |^|2=^t^, (1.2) 

where ^ is one of the scalar multiplets of the GUT-type sector of the system, whose ex- 
pectation value plays the role of an inflaton and which has a strong non-minimal curvature 
coupling with ^ > 1. Here, Mp = mp/VSTT ^ 2.4 x 10^^ GeV is a reduced Planck mass, 
A is a quartic self-coupling of ^, and f is a symmetry breaking scale. 



The motivation for this model was based in the early papers [0, [T^ on the 

observation that the problem of an exceedingly small quartic coupling A ~ 10^^'^, dictated 



by the amplitude of primordial scalar (density) perturbations 14, 15 1, can be solved by 
using a non-minimally coupled inflaton with a large value of ^, because the CMB anisotropy 
AT/T ~ 10^^ is proportional to the ratio \/A/^, rather than to \/A as in the minimal case. 
Therefore, a small AT/T can be obtained even for A close to unity (but still small enough 
to justify perturbative expansion in A) if ~ 10^. 



Later the model (1.1) with the GUT- type sector of matter fields was used to generate 
initial conditions for inflation Q within the concept of the no-boundary |jl^ and tunneling 
cosmological states [17, 18, The quantum evolution with these initial data was 



considered in [11, 21 1. In particular, it was shown that quantum effects are critically 
important for this scenario. 

A similar model but with the Standard Model (SM) Higgs boson <P playing the role of 
an inflaton instead of the abstract GUT setup of |l|, |ll| was suggested in . This work 
advocated the consistency of the corresponding CMB data with WMAP observations. 
However, it disregarded quantum effects^ which were later taken into account in at 
the level of the one-loop approximation. Application of the general method of jl], |ll[ 
to the Standard Model Higgs inflation has led to the lower bound on the Higgs mass 
Mh ^ 230 GeV, originating from the observational restrictions on the CMB spectral index 
m. However, this pessimistic conclusion, which contradicts the widely accepted range 
115 GeV< Mh < 180 GeV, did not take into account 0(1) contributions due to running 
logarithms in the quantum effective action, which can qualitatively improve the situation^. 
This was nearly simultaneously observed in [^, ^, |5| where the renormalization group 
(RG) improvement of the one-loop results of was considered and it was found that the 
Higgs-mass range compatible with the CMB data almost coincides with the conventional 
one. 

In contrast to the numerical and, therefore, not entirely transparent analysis of @, |5|, 
we develop here a method which retains the analytical formalism of Q by incorporating 
the technique of running coupling constants. This supports the results of and, moreover, 
explicitly reveals a mechanism which brings the CMB compatible range of the Higgs mass 
closer to the conventional Standard Model domain. This mechanism also explains the 
perfect efficiency of the perturbation theory at energy scales of the inflationary stage; since 
this is analogous to asymptotic freedom in QCD, we call it here asymptotic freedom, too. 
This mechanism is mediated by the effect of longitudinal (virtual) gravitons; gravitation 
thus serves as a regulator making the Standard Model sector softer. 

The organization of the paper is as follows. In Sects. 2 and 3 we discuss the notion 
of anomalous scaling and recapitulate the results of the one-loop approximation of Q . In 
Sect. 4 we develop the RG improvement for the effective action of the theory, present beta 
functions for running couplings with a special emphasis on loop contributions of Higgs 



^This was caused by the use of a special field-dependent renormalization scheme [g2| suppressing the 
contributions of running logarithms in the effective potential. 

^The possibility of such contributions was, in fact, mentioned in the paper Q, but the sign of their effect 
was not predictable at the moment of publication. 



particles and Goldstone modes in the presence of the non-minimal coupling to graviton 
fn\ , 1^, ^. In Sect. 5 we prove the applicability of the analytic method of |^] for the 
CMB dynamics at the inflationary stage, which is based on the shortness of this stage 
compared to the post-inflationary evolution to the present-day electroweak vacuum. Sect. 
6 contains results of numerical simulations for the RG flow of coupling constants, the 
anomalous scaling A{t) in particular, which impose CMB and Standard Model bounds on 
the Higgs mass and reveal the asymptotic-freedom mechanism based on the behavior of 
A{t) and X{t). In concluding Sect. 7, we discuss the limitations of the obtained results 
associated with gauge and parametrization dependence of the underlying formalism and 
suggest the method of their resolution. We also compare in the context of these issues our 
conclusions with the results of and [|| . We show that the RG treatment of the model is 
essentially more sensitive to the CMB bounds than in |p, especially at the upper limit of the 
Higgs mass range ~ 185 GeV. We also find that the RG improvement makes the scenario 
realistic in which the initial conditions for inflation are generated from quantum cosmology 
in the form of a sharp probability peak of the tunneling wavefunction |17, 18, 19, 20|. 
The Appendix contains the derivation of SM beta functions in the presence of a strong 
non-minimal curvature coupling of the Higgs field. 



2. Anomalous scaling parameter in theories with a large non-minmal cur- 
vature coupling 

The usual understanding of non-renormalizable theories is that renormalization of higher- 
dimensional operators does not effect the renormalizable sector of low-dimensional opera- 
tors, because the former ones are suppressed by powers of a cutoff - the Planck mass Mp 
p^ . Therefore, beta functions of the Standard Model sector are not expected to be modi- 
fied by graviton loops. The situation with the non-minimal coupling is more subtle. Due 
to the mixing of the Higgs scalar field with the longitudinal part of gravity in the kinetic 
term of the Lagrangian ( |1.1| ) , an obvious suppression of pure graviton loops by the effective 
Planck mass, Mp + » Mp, for large ^ proliferates to the sector of the Higgs field, so 
that certain parts of beta functions, which induce Landau poles, are strongly damped by 
large ^ Therefore, running coupling constants like X{t) turn out to be small or remain 
finite at the inflation scale. In particular, a special combination of coupling constants A 
which we call anomalous scaling |l|] becomes very small and reduces the CMB-compatible 
Higgs-mass bound. The importance of A follows from the fact observed in ||l|, 11, ^ that 



due to large ^, quantum effects and their CMB manifestation are universally determined by 
A. The nature of this quantity, which was initially introduced in the context of a generic 
gauge theory non- minimally coupled to gravity Q, is as follows. 

Let the model contain in addition to ( |l.lD also a set of scalar fields Xj vector gauge 
bosons and spinors ip, which have a typical interaction with ^ dictated by the local 
gauge invariance. If we denote by if the inflaton - the only nonzero component of the 
mean value of ^ in the cosmological state, (p"^ = (p^^, then the non-derivative part of the 



interaction of these particles with (/? can be schematically written down as 

i^int = -Yl ^^xxV - Yl - y^^^^^ (2-1) 

X A tp 

where A^^, qa and are the relevant constants of quartic, gauge, and Yukawa couplings. 

On the background of a slowly varying inflaton these terms generate masses m{ip) 
for all these particles, which are proportional to ip, m{(p) ~ (p, and which in turn generate 
at one-loop order the Coleman- Weinberg potential 

(±1) ^-h^ln ^ = :rip^\n-^ + ... . 2.2) 

^ ' 647r2 u2 I287r2 ^ ^ ^ 

particles 

Here, the summation over particles includes the statistics as well as the sum over polariza- 
tions, and the overall coefficient A schematically reads 



2 
A 



E^x + 3E^^A-4E4h (2-3) 

C A ip 



where the coefficients 1, 3 and 4 in front of the sums are the number of degrees of freedom 
of a real scalar field, a massive vector field, and a charged Dirac spinor^. This quantity is, 
of course, well known, because as a coefficient of the logarithm it comprises the ultraviolet 
renormalization of A and constitutes the conformal anomaly, or anomalous scaling, asso- 
ciated with the normalization scale fj, in ( |2.2| ). It includes the one- loop contributions of 
all particles except the graviton and the inflaton fleld (p itself, because their contributions 
are strongly suppressed by inverse powers of » 1 |11, |5| - a property which will be 
discussed in much detail below. 

The role of (|2.3|) is two-fold. First, in the context of the no-boundary and tunnel- 



ing initial conditions 17, [T^, ^ extended to the one-loop level it appears in 
the distribution function of the quasi-de Sitter cosmological instantons p{(p) describing a 
quantum distribution of cosmological models with different initial <p, 

p{^) ~ (^^J . (2.4) 

This scaling (which gave rise to the name "anomalous scaling") makes this distribution 
normalizable at — > oo for positive A |2^ and, moreover, for the case of the tunneling 
cosmological wavefunction it generates a sharp probability peak in p{^) at the following 
value of the infiaton field: the quantum scale of inflation derived in Q, (/j^ = QAir'^Mp/S.A. 

Secondly, the anomalous scaling for ^ ^ 1 determines the quantum rolling force in 
the effective equation of the inflationary dynamics [0, 21 1 and, consequently, yields the 



parameters of the CMB generated during inflation |^ . These parameters explicitly depend 



^For a generic model, (A^, Sai J/V") comprise coupling constant matrices which give rise to mass matrices 
and their traces in the actual expression for A. 



on A and the e-folding number of the first horizon crossing by the primordial cosmolog- 
ical perturbation of a given wavelength. This dependence for the spectral index Ug, which 
belongs to the interval 0.94 < ns{ko) < 0.99 (the combined WMAP+BAO+SN data at the 
pivot point ko = 0.002 Mpc^^ corresponding to ~ 60 |26, gives the range of the 
anomalous scaling — 12 < A < 14 |Q. 

On the other hand, in the Standard Model A is expressed in terms of the masses of 
the heaviest particles - boson, Z boson and top quark, 

2 ^ 2 2 2 ^ r 2 , /2\ 2 2^22 ^oc;^ 

mw = ^9 V , mz = -{g +g ) 93 , = - v? , (2.5) 
and the mass of three Goldstone modes 

ml = YM = x{^''-v')^Xip\ (2.6) 

Here, g and g' are the SU{2) x U{1) gauge couplings, gg is the SU{2>) strong coupling 
and yt is the Yukawa coupling for the top quark. We work in the transversal Landau 
gauge with no contribution to ttoq from the gauge-fixing term. At the inflation stage with 
4P' ^ v^, the Goldstone mass is non- vanishing in contrast to its zero on-shell value in 
the electroweak vacuum [^]. It is again important to emphasize that the Higgs particle 
itself with the mass 



2^ = V"{^) = A(3(/92 - i;2) ~ 3A(/j2 (2.7) 



does not contribute to A, because its contribution is suppressed at the inflation scale by a 
small factor ~ l/C'^; this is shown in the Appendix. 
Equation (|2.2D then gives the expression 

A = A (2/ + [g^ + g'^f - l&yt) + 6A. (2.8) 



In the conventional range of the Higgs mass 115 GeV< M^j < 180 GeV |29(| , this quantity 
is at the electroweak scale in the range —48 < A < —20, which strongly contradicts the 
CMB range given above and violates as well the normalizability of the distribution function 
( |2.4| ). These two ranges can be brought together only by the price of raising A ~ Mfj/2v'^, 
that is, by increasing the Higgs mass to 230 GeV Q. 

Salvation comes, however, from the observation that the RG running of coupling con- 
stants is strong enough and drives A to a range compatible with the CMB data and the 
conventional range of the Higgs mass. Thus, the formalism of |^] stays applicable but with 
the electroweak A replaced by the running A(t), 



A{t) = ^(2/(t) + {9'{t)+g'\t)y - 16yfit)) +6A(t). (2.9) 
Here, t = ln{(p/fj,) is the running scale of the RG improvement of the effective potential 



1 3C ] , where /x is a normalization point which we choose to coincide with the top quark mass 
II = Mt (we denote physical (pole) masses by capital letters in contrast to the running 
masses (|2.5|) above). In the following we shall present the details of this mechanism. 



3. One-loop approximation 



The inflationary model with a non-minimally coupled Higgs-inflaton was considered in the 
one-loop approximation in Q. The low-derivative part of its effective action (appropriate 
for the inflationary slow-roll scenario), 



<fxg^'^-V{^) + U{^)R{g^,) 



contains coefficient functions which in the one-loop approximation read as 



1287r2 

.2 



Cln^ + D 



G(^) = 1 + 



1 



327r^ 



(3.1) 

(3.2) 
(3.3) 
(3.4) 



These functions contain numerical coefficients A,C,D,F,E determined by contributions 
of quantum loops of all particles and include the dependence on the UV normalization 
scale /i. For ^ 1 the inflationary stage and the corresponding CMB parameters critically 
depend only on the anomalous scaling A defined by ( |2.2| ) ||l|, 11, |3| and the part of C 
which is linear in ^. As shown in Appendix A, the latter is contributed by the Goldstone 
mass niQ and equals 

C = 3^A + 0(C°). (3.5) 

In fact, this expression yields a quantum correction to the non-minimal coupling as can 
be seen from (3.3). Other coefficients, the normalization scale fx inclusive, are irrelevant in 
the leading order of the slow roll expansion. 

Inflation and its CMB are easy to analyze in the Einstein frame of fields g^^, (f, which 
is related to the Jordan frame used in (p.ip by the equations 



2^((p) 

M2 



?72 



(3.6) 



The action (3J.) in the Einstein frame, S[gfj_y^0\ = S[gfj_y^(p\^ has a minimal coupling, 
tj = Mp/2, a canonically normalized inflaton field, (5 = 1, and a new inflaton potential. 



V{0) = 



Ml 



(3.7) 



In view of (jSJ), (|3J) and (|3 
^ » 1, this potential reads 



V 



AM|, 



at the inflation scale with > Mpl\/\ ^ v and 

(3.8) 



+ — ^ In 




where the parameter Aj represents the anomalous scahng (^^) modified by quantum 
corrections to U ~ the renormahzation of the non-minimal term (^]^) , 

Ai = A-12X = ^ (2/ + {g^ + ' - 16^^^) - 6A. (3.9) 

This quantity - which we will call inflationary anomalous scaling - enters the expressions 
for the standard slow-roll parameters, 

^ A / Ml Ai\^ , , 

MlfV__iMl 

and ultimately determines all the characteristics of inflation. In particular, the smallness 
of e yields the range of the inflationary stage 99 > (/9cndi terminating at a value of e which 
for elegance of the formalism we chose to be Send = 3/4. Then the inflaton value at the exit 
from infiation equals (^cnd — IMpj under the natural assumption that the perturbation 
loop expansion is applicable for A//647r^ <C 1. The duration of inflation which starts at 99 
takes in units of the scale factor e- folding number N a particularly simple form |Q , 

^ = e-- 1, (3.12) 

"Pi 

V>i = , (3-13) 
c.Ai 

where the quantum scale of inflation ipj and a special parameter x containing the e-folding 
number, 

x = ^ (314) 

directly involve the anomalous scaling A/ which is essentially a quantum quantity. 

This relation determines the Fourier power spectrum for the scalar metric perturbation 
C, A^(A;) =< k'^Ck V /24:Tt'^ Mpi, where the right-hand side is taken at the first horizon 
crossing, k = aH, relating the comoving perturbation wavelength k"^ to the e-folding 
number N, 

The CMB spectral index = 1 + d In A^/d In A; = 1 — 6i + 2fj, the tensor to scalar ratio 
r = 16i and the spectral index running a = dUg/dlnk ~ —dus/dN correspondingly read 
as 

2 X , , 



Note that for |x| <C 1 these predictions exactly coincide with those ^] of the f{R) = 
Mp(i? + i?^/6M^)/2 inflationary model with the scalar particle (scalaron) mass M = 

MpVx/V^i- 

With the spectral index constraint 0.94 < ns{kQ) < 0.99 at the pivot point kQ = 0.002 
Mpc^^ corresponding to ~ 60 |2(:, 27] these relations immediately give the range of 
anomalous scaling — 12 < A/ < 14. As mentioned above, this contradicts the Standard 
Model bound —48 < Aj < —20 existing in the conventional range of the Higgs mass 115 
GeV< Mh < 180 GeV 

The limitation of this approximation follows from the fact that, despite the smallness of 
Ai /64Tr^, the running of the logarithm from the electroweak scale (p = v (where the Higgs 
and other particle masses are determined) to the inflation scale (/9 ~ \ipj\ is big enough and 
changes the predictions by 0(1) factors. Indeed, with the values of ~ lO'' and \Az\ ~ 10 
dictated by the results of Q, the estimate for the logarithmic running ln(|99||/f^) ~ 60 
yields the contribution 

ln^~2 = 0(l), (3.19) 



647r2 



which requires resummation within the RG improvement. The necessity of this improve- 
ment was mentioned in Q, and it was realized in @, |5|. Below we develop a combined 
numerical and analytical scheme of this improvement, which recovers the above formalism 
of Eqs. (3.12|)-(3.1^) incorporating running coupling constants. 



4. RG improvement 



According to the Coleman- Weinberg technique [30|, the one- loop RG improved effective 
action has the form (13. 11) with the coefficients 



uiip) = ^(^Mi,+mz\t)ip^), 

G{^) = Z\t). 



(4.1) 

(4.2) 
(4.3) 



Here, t = ln((/?/Aft) is the running RG scale^, and the running couplings \{t), ^{t) and the 
field renormalization Z(t) incorporate a summation of powers of logarithms and belong to 
the solution of the RG equations 



dm 

dt 
dZ 

'dt 



■ Pg^^ gi = {\L9,9',9s,yt), 

iz 



(4.4) 
(4.5) 



^These estimates were obtained in for the parameter A disregarding the contribution of Goldstone 
modes (+6A and — 6A terms in (2.6) and (3.9), respectively), but numerically with A < 1 the inclusion of 
this contribution leads to qualitatively the same conclusions. 

^Application of the Coleman-Weinberg technique removes the ambiguity in the choice of the RG scale 
in cosmology - an issue discussed in |3]{|. 



for the full set of coupling constants in the "heavy" sector of the model with the correspond- 
ing beta functions and the anomalous dimension 7 of the Higgs field. ^ This sector of 
the Standard Model includes the SU{2) x U{1) gauge couplings g and g', the SU{3) strong 
coupling gs, and the Yukawa coupling yt for the top quark. 

An important subtlety with these /3 functions is the effect of the non-minimal curvature 
coupling of the Higgs field. For large ^ the kinetic term of the tree- level action has a 
strong mixing between the graviton and the quantum part of the Higgs field a on the 
background if. Symbolically it has the structure 



which yields a propagator whose elements are suppressed by a small 1 /^-factor in all blocks 
of the 2x2 graviton-Higgs sector. For large ip ^ Mp/\/^, the suppression of pure graviton 
loops is, of course, obvious because the effective Planck-mass squared exceeds by far the 
standard Planck-mass squared, Mp + ^(/j^ ^> Mp. Due to the mixing, this suppression 
proliferates to the full graviton-Higgs sector of the theory. To make this statement quan- 
titative, one can go over to the Einstein frame (|3.6| ) in which the propagator is diagonal 
and canonically normalized in the space of perturbations of g^^j and ip (in the background 
covariant DeWitt gauge) . The above propagator can be uplifted back to the Einstein frame 
by matrix multiplication with the Jacobian matrices d{g^i,,ip)/d{ga/3,^)- Its Higgs-Higgs 
block, in particular, gets suppressed by the factor s{ip) 



00 J g^'^in-ml) g^/^n-ml)' 

^^'^^ - Gc/ + 3?7'2 M2 + (ee + i)e</'2 ' ^^■') 

where we took into account the scaling g^^'^{\I\ — rh'jj) = {2U/Mp)g^^'^(n\ — m'jj) under 



the conformal transformation (3.6) (disregarding spacetime gradients of (p and bearing in 
mind conformal rescaling of particle masses j^). The origin of this suppression factor s{ip) 
within the original Jordan frame is clearly shown in Appendix A. 

This mechanism was first understood in P] when estimating the CMB generation in a 



non- minimal inflationary model. It was considered in |35| in the context of the generalized 



RG approach to nonrenormalizable theories. Also it justifies the omission of graviton 



loops [l^, ^ and modifies the beta functions of the SM sector of this theory at the 



high energy scale relevant for infiation. This modification is due to the suppression of 



every Higgs propagator (4^) by the factor s{ip) |5|, which is very small for ip » Mp/y/^, 



s ~ 1/6^, but tends to one in the vicinity of the electroweak scale v <C Mp/^. Such a 



®Note that running couplings in the effective potential theory are the characteristics of the RG equation 
in partial derivatives for this potential, rather than the usual running coupling constants directly run by /3's. 
Thus the right-hand sides of RG eqs.([4.^)-(4.?: ) should contain the factor 1/(1 — 7) which we discard 
in the one-loop RG improvement. 



modification justifies, in fact, the extension beyond the scale Mpj^ interpreted in |36, 37 1 
as a natural validity cutoff of the theory'^. 



Of course, the factor s{}f) makes the counterterms of the theory nonpolynomial in 
y?, and the theory as a whole becomes nonrenormalizable. One might think that this 
completely invalidates the construction of the usual RG improvement and requires the 
generalization of the latter to RG with an infinite set of charges (perhaps of the functional 
nature like in Q). However, there exists a shortcut to a simpler formulation based 
on the structure of the factor (|4.7| ). In the logarithmic scale of the variable t = ln{ip/fj,), 
strongly compressing the transition domain Mp/£^ ^ f ^ ^p/V^ for •s(v) from 1 to 
1/6^ ~ 0, this factor looks very much like a step function 9{tQ — t), to ~ ln(<y9o//i) = 
So it is nearly constant for (f < (fo and (p > (fo- Therefore, in both 
phases of such an approximation the theory is renormalizable but has different running of 
coupling constants. This qualitatively justifies the use of beta functions modified by s- 
factors smoothly interpolating between these two phases. Below we will see that numerical 
results for smooth and step function s-factors nearly coincide with one another. 

There is an important subtlety with the modification of beta functions, which was 
disregarded in |||] (and the first version of this paper). The mnemonic rule of associating 
the factor s{ip) with every propagator of the Higgs multiplet, as suggested in turned 
out to be incorrect. Goldstone modes, in contrast to the Higgs particle, are not coupled to 
curvature, and they do not have a kinetic term mixing with gravitons [^. Therefore, their 
contribution is not suppressed by the s-factor of the above type. Separation of Goldstone 
contributions from the Higgs contribution leads to the following modification of the one- 
loop beta functions, which is essentially different from that of Q (cf. also [^]): 



Px-- 



X 



167r2 

JL 
167r2 

■ 

' 167r2 



{l8s'^X + A{t)) -47A, 
(l + s')A-27e, 



J2 



(4.8) 
(4.9) 
(4.10) 



^The notion of a cutoff depends on the type of perturbation theory used. Typically it is determined 
for the case when all dimensional quantities - fields and their derivatives - are considered on equal footing 
and treated perturbatively. However, a nonperturbative summation of powers of the field for low-derivative 
processes (achieved by the conformal transformation to the Einstein frame) allows one to go beyond such 
a cutoff, and this is exactly the situation with the scales above A/p/^, see e.g. jssj. 

*This is analogous to the approach of Q where the model was approximated by matching the usual 
low-energy SM phase with the chiral phase of the SM at the infiation energy scale. This picture, however, 
takes place in the Einstein frame of the theory and is aggravated by the problem of transition between two 
different parameterizations of one quantum theory — the Cartesian coordinates in the space of the Higgs 
multiplet in the low-energy phase versus spherical coordinates in the chiral phase of the SM. This transition 
is physically nontrivial and is discussed below in more detail. 



and 

B--^1^J1- (4 11) 

= (4.12) 

Here, the anomalous dimension of the Higgs field 7 is given by a standard expression in 
the Landau gauge, 

the anomalous scaling A{t) is defined by (|2.9| ), and we have retained only the leading 
terms in ^ ^ 1. It will be important in what follows that this anomalous scaling contains 
the Goldstone contribution 6A, so that the full Px in ([4.8|) has a A^-term unsuppressed 
by s{ip) at large scale t = ln{ip/fi). In Appendix A we derive the Goldstone and Higgs 
contributions to f3x and which demonstrates their different suppression mechanisms 
under non-minimal coupling with curvature. 

5. Inflationary stage versus post-inflationary running 

As we will see, the infiationary stage in units of Higgs-field e-foldings is very short. This 
allows us to make one more approximation: we shall consider the solutions of the RG 
equations only up to terms linear in At = t — tend = ^^{v/Vend)^ where we choose as the 
initial data point the end of infiation tend- This approximation will be justified later in 
most of the Higgs-mass range compatible with the CMB data. 



Thus, we use beta functions (4.8) and (|4.9|) with s = to obtain 



A(t) = Aend (1 - 47endAt + At) , (5.1) 

m = ^end (l - 27cndAt + ^^At) . (5.2) 

Here, Aend; Tend) '^end are determined at tend and A^nd = ^(*end) is also the particular value 
of the running anomalous scaling ( |2.9| ) at the end of inflation. 

On the other hand, the RG improvement of the effective action ([4.1|)-( [43| ) implies 
that this action coincides with the tree-level action with running couplings as functions of 
t = ln{ip/ fi) for a new field 

(p = Z{t)ip. (5.3) 

This is because the running of Z(t) is slow, 5^0 = Z(l + 'j)d^(p ~ Zd^ip, and the kinetic 
term of the effective action in terms of (j) gets canonically normalized. Then, in view of 



(4J)-(4^) the RG improved potential for this field takes at the inflation stage the form 



,> iKVv .,4 /, 2Ml AiiU^i, ^ . 



where for the same reason we disregarded the running of Z in ln(Z/Zg^(i). This is nothing 



but the one-loop potential (3.8) for the field (j) with a particular choice of the normalization 
point /i = (j)end and the couplings replaced by the values of the running ones at tend- 

This means that the formalism of |^ can directly be applied to determine the pa- 
rameters of the CMB. They are mainly determined by the anomalous scaling A/, but 
now this quantity should be taken at tend rather than at the electroweak scale t = 0. 
This can solve the problem of matching the electroweak range of the anomalous scaling 
-48 < Ai{0) < -20 (corresponding to 115 GeV < Mh < 180 GeV) with the range re- 
quired by the CMB data, -12.4 < Ai{tend) < 14.1 [|], because A/(tend) 7^ Ai{0) due to 
the running. 

To find A/ (tend) and other couplings at tend we have to obtain the RG flow interpolating 
between the t = and tend — 33. In contrast to the inflationary stage, the post-inflationary 
running is very large, because of its long duration, and requires numerical simulation. In 
order to specify more precisely the initial conditions at the top-quark scale, t = 0, we 
take into account the pole mass matching scheme relating the observable physical masses 



Mh and Mt to the running masses mn = '\/2X{t)v and (2^) (and relevant couplings 



1 41, 43 1 . One might think that these subtleties essentially exceed the precision of 
possible SM implications in early cosmology. However, the instability bound on the Higgs 
mass and its lower bound from the CMB data, which turn out to be very close to one 
another, are very sensitive to the pole mass matching effect, so we include them in numerical 
simulations. 

To have a better comparison with ||5|, Q, we fix the t = initial conditions for the 
RG equation (|4.4| )-( p^ at the top-quark scale Mt = 171 GeV. For the weak interaction 



constants g,g' and the strong interaction constant gs, they read 

5^(0) = 0.4202, g'2(0) = 0.1291, (/^(o) = 1.3460, (5.5) 

where 5^(0) and ^'^(O) are obtained by a simple one- loop RG flow from the conventional 
values of a{Mz) = g^/^7r = 0.0338, a'{Mz) = ^'V^tt = 0.0102 at the M^-scale, and 
the value g1{Q) at Mt is generated by the numerical program of |44]^. For the Higgs 
self-interaction constant A and for the Yukawa top quark interaction constant yt the initial 



conditions are determined by the pole mass matching scheme originally developed in |41, 42 1 



and presented in the following form in the Appendix of [43|: 



j,f2 

A(0) = -|(1 + 2Ah(Mh)), (5.6) 
yt(0) = ^^(l + Ai(MH)). (5.7) 

V 

Here, the electroweak vacuum expectation value for the Higgs field is f = 246.22 GeV, and 
/S.h{Mh) and A(M//) comprise relevant mass operator corrections in the effective Higgs 
and top quark propagators (cf. the relation (|2.5D for masses without these corrections), 
which together with Mh depend on Mt, Mz and the Weinberg angle. 



The analytical algorithm of transition between diflterent scales for gl was presented in [M5| 



The initial condition ^(0) is not directly known. It should be determined from the CMB 
normalization condition for the amplitude of the power spectrum (|3.15| ) ~ 2.5 x 10"^ 
at the pivot point ko = 0.002 Mpc"^ |2|, |2^ which we choose to correspond to ~ 60. As 



quantum perturbations of the observable ip and the new field ( ^.31 ) are obviously related by 
Ci/p = C(f,/Z, Eq. ( 3.15| ) immediately yields the following estimate on the ratio of coupling 



constants: 

-2--^ ~ 0.5 X 10 9 5.^ 

at the moment of the first horizon crossing for N = 60, which we call the "beginning" of 
inflation and label by ti^.^^ This moment, in turn, can be determined from the equations 



( 3.12 )-( 3.14 ) and ( 3.13 ) relating the value of the inflaton field at this moment (pi^i to the 
e-folding number and anomalous scaling A/(tin). With ipin = Mjexp(tin), this relation 
takes the form , ^ , , ^ 

im = ln— + -ln— + -ln . 5.9) 

Mt 2 34in 2 Xin 



The RG equations ( |4.4| ) for the six couplings {g, g' , gs,yt, ^,0 with five initial condi- 
tions (|5.5D - (|5.7| ) and the final condition (|5.8D at tin defined by Eq. (|5.9|) uniquely determine 
the RG flow for given values of the Higgs and top quark masses. Boundary conditions for 
this flow are rather involved, because they go beyond a usual Cauchy problem and even the 
deflnition of the final moment ( |5.9D includes the whole history from t = to t = tin, that 
is, they are nonlocal in t. By iterations, however, this problem can be numerically solved 
with Mathematica; essential simplifications originate from the fa-ct that Xjn turns out to be 
0(1) (more precisely ~ 2, see below) and that the last term in ( |5.9D almost always can be 
discarded from the expression for tin on top of the contribution ln{Mp / Mt) ~ 37. 

The RG fiow covers also the inflationary stage from the chronological end of inflation 
^end to tin (remember that the arrows of the RG parameter t and the physical time are 
opposite). As was mentioned above, at the end of inflation we choose the value i = 



3/4 of the slow roll parameter ( ^.10 ), and (fend = ^P\l 4/3^cnd under the assumption 
A7(tcnd)/647r^ ^ 1 which will be justifled below. Therefore, 

Mp 1 4 

W = ln-rf + :.ln-— . (5.10) 

Thus the duration of inflation in units of inflaton field e-foldings tin — ^end = ^'^{^m/ ^cnd) 
is very short relative to the post-inflationary evolution tend ~ 35, 

tin - tend = \\n N + \\n + \\n !^P^HL^ 1 In iV ~ 2, (5.11) 

2 Send 2 ■^in ^ 

where we took into account that ^in — ■^cnd a.nd x = 0(1) . This is, of course, typical for 
large-field inflationary models with ~ ip^ /Mp. 

This estimate gives the range of our approximation linear in logarithms. The lineariza- 



tion, say in (5.1), implies the bound | A7(tend)|^t/167r <C 1, where At < tin— tend — lnA^/2, 



^"No modification of tlie spectral index due to Z occurs in the one-loop RG running, because tfie difference 
t -nt = 2'Ydt/dN = (7y4//487r2)e"=/(e'= - 1) belongs to the two-loop order. 



and this approximation holds for 



\Ai{t. 



end J 



167r2 



InA^ 



0.5. 



(5.12) 



6. Numerical analysis 

The running of A{t) strongly depends on the behavior of A(t). It is well known that for small 
Higgs masses the usual RG flow in SM leads to an instability of the electroweak vacuum 
caused by negative values of X{t) in a certain range of t, see e.g. |46, and references 
therein. The same happens with the modified RG for the non-minimally coupled Higgs 
field considered here. 
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Figure 1: Running X{t) for five values of the Higgs mass above the instability threshold. Dashed 
curves mark the boundaries of the inflation domain icnd < i < ^in- 
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The numerical solution for X{t) in the setting of the previous section is shown in Fig.l 
for five values of the Higgs mass and the value of top-quark mass Mt = 171 GeV. The 
lowest one corresponds to the critical value 

Mjj'"* ~ 134.27 GeV. (6.1) 

This is the boundary of the instability window for which A(t) bounces back to positive values 
after vanishing at tinst ~ 41.6 or (fmst ~ 80Mp}^ It turns out that the corresponding ^(t) 
is nearly constant and is about 5000 (see below), so that the factor {\i.7\) at tinst is very 
small, s ~ 1/6.^ ~ 0.00005. Thus the situation is different from the usual Standard Model 
with s = 1, and numerically the critical value turns out to be higher than the known SM 
stability bound ~ 125 GeV ||. 

^^An over-Planckian scale of (p does not signify a breakdown of the semiclassical expansion, because the 
energy density ~ 10~^^ Mp stays much below the Planckian value, see below. 



-14- 



Fig.l shows that near the instabihty threshold Mr = M]^''* the running couphng A(t) 
stays very smah for all scales t relevant to the observable CMB, ( |5.9D -( ^?10|) . This follows 
from the fact that the positive running of \{t) caused by the term (18s^ + 6)A^ in /3a, (0), 
is much slower for s ^ 1 than that of the usual SM driven by the term 24A^. For larger 
Mh this suppression in the term ISs^A^ is responsible for the origin of the bump on the 
plots at t ~ 26 where the value of s drops from 1 to and A(t) continues growing but with 
a slower rate. 

60 ' ' ' \ 
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Figure 2: Running anomalous scaling for the critical Higgs mass (the red curve with a vertical 
segment at the singularity with tinst ^ 41.6) and for two masses in the stability domain (blue and 
green curves). 



The smallness of A improves, of course, the efficiency of the perturbation theory. On 
the other hand, at the critical value of the Higgs mass we encounter a singularity of both 
anomalous scalings A{t) and Ai{t) at the instability point tmst where A(tinst) = 0. This 
is depicted in Fig. 2 for Ai{t) by the red curve corresponding to the slightly overcritical 
value 134.271 GeV (vertical red line and quasi-hyperbolic curves to the left and to the right 
of it). Other overcritical curves with small Mh > Mf^^^ run through zero and change the 
sign of Ai{t) from negative to positive. For larger Mh, when A also becomes larger, the 
term — 6A in Eq. ( ^^ ) for A/ does not let it get positive at any t, and Aj starts decreasing 
after reaching some maximal negative value (like for the green curve of Mh = 180 GeV). 

An important observation is that for all Higgs masses in the range M]^*^* = 134.27 GeV 
< Mh < 185 GeV the inflation range tend < t < tin is always below the instability value 
iinst = 41.6 (numerics gives that tin < 34.8), so that from Fig. 2 A/(t) is always negative 
during inflation. Its running depicted in Fig. 2 explains the main difference from the results 
of the one-loop calculations in Aj{t) runs from big negative values A/(0) < —20 at the 
electroweak scale to small but also negative values at the inflation scale below tinst- This 
makes the CMB data compatible with the generally accepted Higgs mass range. Indeed, 



the knowledge of the RG flow immediately allows one to obtain A/ (tend) and Xend and 
thus find the parameters of the CMB power spectrum ( 3.16| )-( |3.18| ) as functions of Mh- 
The parameter of primary interest - the spectral index - is given by Eq. with 
X = a^end = ^A/(^end)/487r^ and depicted in Fig. ^ Even for low values of the Higgs 
mass above the stability bound, falls into the range admissible by the CMB constraint 
existing now at the 2a confidence level (based on the combined WMAP+BAO+SN data 

Hi), 

0.94 < n,(A;o) < 0.99. (6.2) 
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Figure 3: The spectral index Ug as a function of the Higgs mass Mr for three values of the top 
quark mass. 

The spectral index becomes too small (that is, dropping below 0.94) only for large Xend 
or large negative ^/(tend)) which happens only when Mh either approaches the instability 
bound or exceeds 180 GeV at the decreasing branch of the Ug graph. Thus, we get lower 
and upper bounds on the Higgs mass, which both follow from a comparison of the lower 
bound of the spectral index in ( |6.2D with (3.16). Numerical analysis for the corresponding 
Xend — —1.4 gives for Mt = 111 GeV the following range for a CMB-compatible Higgs 
mass: 

135.62 GeV < Mh < 184.49 GeV. (6.3) 



Both bounds belong to the nonlinear domain of the equation ( ^.16 ) because their 
relevant Xend = —1-4 < —1. However, their calculation is still in the domain of our 
linear in logs approximation, because the quantity A/(tend)/167r^ ~ —0.07 satisfies the 
restriction ( 5.12 ) for Xend = —1-4. For a smaller |A7(tend)| the expression ( |3.16| ) for 
be linearized in x and takes a particularly simple form. 



n, can 



, 2^ A/ (tend) 

~ N 



A/ (tend) « 



487r^ 



(6.4) 



It is applicable in a wide range of the Higgs mass in the interior of the domain ( |6.3D . 

As we see, the upper bound on in (|6.2| ) does not generate any restrictions on Mh, and 
it will not effect the M//-range unless it will be lowered down by future CMB observations 
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Figure 4: The Einstein frame effective potential for the instabihty threshold Mf^^'^ = 134.27 GeV. 
A false vacuum occurs at the instability scale t-mst — 41.6, (p ^ 80Afp, which is much higher than 
the inflation scale < (^in — 0.04Mp. A hypothetical inflation domain (ruled out by the lower Ug 
CMB bound and the requirement of the positive slope of ^ at i < tin) is marked by dashed lines. 



to about 0.964 — the top of the n^-graphs in Fig. ^. The lower CMB bound in (6.3) is 
slightly higher than the instability bound M]^^* = 134.27 GeV. This bound depends on 
the initial data for weak and strong couplings and, what is even more important, on the 
top-quark mass Mt which is known with less precision. The bound M]^^* given above was 
obtained for Mt = 171 GeV. Below we consider basically this value of the top-quark mass. 
Results for the neighboring values Mt = 171 it 2 GeV are presented in Fig. ^ only to show 
how strongly the plot gets shifted along the Mh axis. The general pattern of this shift 
follows the dependence of M]^*^* on Mt - the instability bound is larger for a larger top 
mass, which can be explained by a negative contribution 
/3a, (E 



4 

-Vt ~ 



-Mt to A, and via A to 



Viability of the inflation scenario also implies that the system should safely evolve from 
the end of inflation to the EW vacuum through the periods of thermalization (reheating), 
radiation and matter domination. The main requirement for this is the possibility to roll 
down from ife^d = Mt exp(tend) to = t> or the positivity of the effective potential slope Q 
(we disregard the tunneling scenario). The shape of the graph of this potential depicted in 
Fig. ^ for the instability threshold confirms the danger of having a negative slope for Mh 
close to M)^*** - the formation of a false vacuum at the instability scale. At higher Higgs 
masses till about 160 GeV we first get a family of metastable vacua at scales > tjnst- An 
example is the plot for the lower CMB bound Mjj = 135.62 GeV depicted in Fig. ^. For 
even larger Mh these metastable vacua get replaced by a negative slope of the potential 
which interminably decreases to zero at large t (at least within the perturbation theory 
range of the model), see Fig. However, these metastable vacua and the negative slope 
of the potential are not dangerous for the inflationary scenario whose scales tin < 34.8 are 
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Figure 5: Inflaton potential at the lower CMB compatible value of Mh- A metastable vacuum 
exists at t ~ 42 which is much higher than the inflation domain at the positive slope of the potential. 
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Figure 6: The succession of effective potential graphs for AP^^^ < Mh < 184.3 GeV showing 
the occurrence of a metastable vacuum followed for high Mh by the formation of a negative slope 
branch. Local peaks of V situated at i = 34 35 grow with Mh for Mh < 160 GeV and start 
decreasing for larger Mh- 



much lower than the metastability scale and belong to the positive slope of V to the left of 
its peak at t = 34 35. Thus, even if the metastable vacuum occurs for low Mh, it exists 
before the inflation stage probed by current CMB observations^^. 

Let us finally focus on the running of ^(t) depicted for five values of the Higgs mass 

^^The existence of this vacuum can perhaps be probed by wavelengths longer than that of a pivotal 
~ 60, but this requires a deeper analysis. 
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Figure 7: Plots of running ^(i). 

in Fig. 1^ starting with the lower bound of the range ( |7.1| ) below. It is very slow for low 
values of the Higgs mass near the instability threshold, which of course follows from the 
smallness of the running A(t) in this domain. Another property of the ^-behavior is that 
the normalization of the power spectrum ( |5.8| ) leads to a value ~ 5000 for small Higgs 
masses, which is smaller than the old estimate ~ 10^ [@, |9|, ||, |ll|, §, |3|. This is caused by 
a decrease of X{t) which at becomes much smaller than A(0) - an observation made in 
Q . This relaxes the situation with the "unnaturalness" problem in this model and brings ^ 
closer to a very subjective borderline between "natural" and "unnatural" values of coupling 
constants. For large Mh close to the upper bound of the CMB range, ^(t) grows like X{t) 
to the Landau pole and eventually goes out of the perturbation theory domain. However, 
it grows much faster than A(t), so that V{ip) ~ A(t)/^^(t) tends to zero, as shown in Fig. 
|6| by the negative slope branches of the potential. 

A final comment of this section concerns the change of numerical results for a "two- 
phase renormalizable" model of the system when the s-function is replaced by the step 
function, discussed in Sect. 4 as a justification for the modified RG scheme. This replace- 
ment leads to a negligible change in numerics - the instability threshold remains nearly the 
same, lower CMB bound on Mh grows by one in the second digit after the decimal point, 
while the upper bound decreases by 1 GeV. 

7. Conclusions and discussion 

We have found that the considered model looks remarkably consistent with CMB observa- 
tions in the Higgs-mass range 



135.6 GeV < Mh < 184.5 GeV, 



(7.1) 



which is very close to the widely accepted range dictated by electroweak vacuum stability 
and perturbation theory bounds. 

Both bounds follow from the lower WMAP bound for the CMB spectral index ns(A;o) > 
0.94 (the combined WMAP+BAO+SN data at the pivot point /cq = 0.002 Mpc"^ corre- 



sponding to ~ 60 e- folds [26, 2^). The lower bound on Mh is very close to the instability 
threshold which in this model equals Mfj^^ — 134.27 GeV. This is higher than the conven- 
tional SM estimate ~ 125 GeV, which is caused by a strong non-minimal back reaction 
of gravity and/or the simplicity of our approximation disregarding higher loop orders. 
Remarkably, the upper bound in (7.1) does not follow from a demand of a valid perturba- 
tion theory, but rather is imposed by the CMB data. This CMB mechanism works in a 
trustable perturbation regime with A(tin) < 2. 

The current upper CMB bound Ug < 0.99 does not impose restrictions on the Higgs 
mass, but it will generate bounds if the observed Ug essentially drops below the "classi- 
cal" value Us = 1 — 2/N ~ 0.967. Conversely, this model will be falsified if the lower 
observational bound on exceeds 0.964 — the top of the n^-graphs in Fig. ^. 

Our approach represents the RG improvement of our analytical results in |^. In fact, 
here we completely recover the analytic formalism of Q for all inflation parameters, which 
only gets modified by the RG mapping between the coupling constants at the EW scale 
and those at the scale of inflation. A peculiarity of this formalism is that for large ^ S> 1 
the effect of the SM phenomenology on inflation is universally encoded in one quantity - 
the anomalous scaling Aj. It was earlier suggested in Q for a generic gauge theory (^), 
and in the SM it is dominated by the contributions of heavy particles - {W^, Z)-bosons, 
top quark and Goldstone modes. This quantity is forced to run in view of RG resummation 
of leading logarithms, and this running raises a large negative EW value of A/ to a small 
negative value at the inflation scale. Ultimately this leads to the admissible range of Higgs 
masses ( |7.1| ) very close to the conventional SM range. 

In fact, this mechanism can be regarded as a kind of asymptotic freedom, for Ai/GAir'^ 
determines the strength of quantum corrections in inflationary dynamics |jll|, ^. Usually, 
asymptotic freedom is associated with the asymptotic decrease of X{t) to zero. Here, this 
phenomenon is trickier because it occurs in the interior of the range ( |7.1| ) and fails near 
its lower and upper boundaries. For small Higgs masses in the range (7.1), as it follows 
from the plots of Fig. ^ and Fig. ^, running A(t) gets very small when approaching the 
inflationary region t ~ 32 35 and generates a large Ai{t) because of its 1/A part induced 
by vector boson and top quark particles. For large masses the coupling A grows towards the 
Landau pole (though much slower than in nongravitating SM because of the suppression of 
the major part of the A^-term of its beta function). Therefore, A/(t) becomes big negative 
again, this time due to the negative contribution — 6A of the Goldstone modes, and falls 
out of the CMB range at Mh ^ 185 GeV. 

Quantum effects are small only in the middle part of ( |7.lD with a moderately small A 
where Ug is close to the "classical" limit 1 — 2/A^ ~ 0.967 for x = NA/ASir'^ <C 1. Here 



It is well known that despite the perturbative range of coupling constants the two- loop RG improvement 
essentially lowers down the EW instability threshold compared to the one- loop RG running |43[|. 



the original claim of Q on smallness of quantum corrections is right, but this smallness, 
wherever it takes place, is achieved via a RG summation of big leading logarithms. 

Qualitatively our main conclusions are close to those of ||5| and |p, though the RG 
treatment in these works is very different from ours. In contrast to them, however, we did 
not try to suggest error bars for the Higgs mass range ( [7. ID and specify corrections due to 
the uncertainty in the top quark mass Mt- We only presented all our results for the middle 
value of the top mass Mt = 171 GeV and demonstrated the sensitivity of the results to this 
choice by two plots of Ug with two neighboring values 169 GeV and 173 GeV (see Fig. |3|). 
These error bars as well as the inclusion of the two-loop approximation are, in our opinion, 
not reliable. Generally speaking, all the conclusions of |^, § and our work still sound too 
good to be true in view of the conceptual problems that have not yet been resolved. These 
problems were partly mentioned in the criticism of Q by Bezrukov and Shaposhnikov Q, 
and now we would like to briefly dwell on them. 

As mentioned in Q , the RG flow of Q depends (via the anomalous dimension of the 
Higgs field 7) on the choice of gauge for local SU (2) x U{1) invariance in the SM. This 
inalienable feature of any RG scheme in SM theory was bypassed in as follows. The 
RG flow of all gauge- independent charges of the SM was used at the post-inflationary 
stage, while at the inflation stage the effect of running renormalization of the Higgs fleld 
(causing gauge dependence of the results) was disregarded due to the flatness of the effective 
potential ~ 1/^. In fact, this argument is misleading, because a small gauge- dependent 
slope of the inflation potential ~ (p non-analytically enters the expression for the e-folding 
number N = J dipH{ip)/ip ~ one of the main ingredients of the inflation parameters C, and 
Ug. Therefore, anyway, this brings an order 0(1) gauge dependence into the expressions 
for these parameters. 

In fact, this is a manifestation of a general problem of the off shell extension in quan- 
tum fleld theory. Unfortunately, we cannot formulate the calculation of CMB parameters 
as on-shell amplitudes or expectation values of physical observables uniquely deflned at 
the quantum level. The power spectrum represents the quantum correlation function cal- 
culated at a special fleld-dependent moment of time - the horizon crossing. By using 
Heisenberg equations of motion (both for linear quantum perturbations and the nonlinear 
background) this function expresses as = t^/247r^Mpf entirely in terms of the inflaton 
potential V{ip) and its gradients [||, ||, H, cf. Eq.(|l^) for e. So the main inflation 
observable is the inflaton potential itself, but its quantum version V'^^ and its off-shell 
extension with {y'^^{ip))' 7^ - the effective potential - is gauge-dependent^^. Moreover, 
off-shell extension of the effective potential (and more generally effective action) also de- 
pends on parametrization of quantum flelds. A particular example of this dependence is 
given by different results for CMB characteristics obtained in different conformal frames 
- the original Jordan frame versus the Einstein one |6|. Another important example 
is the SM renormalization with different parameterizations of the Higgs multiplet. In the 
conventional quantization in Cartesian coordinates for the complex SU(2) doublet or real 
0(4) multiplet the theory is renormalizable and has one set of quantum corrections. In the 

^''in fact, this is a problem of consistent transition from {0\V{(p)\) to 0\y'^^ {{ip))] which brings into game 
the problem of gauge dependence in the formalism of the mean field (<^) and its effective action. 



spherical coordinate system (used, for example, for the description of the chiral phase of 
SM E^, |6|) the theory is non-renormalizable and has another set of quantum corrections 



1 48]. The S-matrices in both parameterizations coincide, but their off-shell effective actions 
are different. 

Both off-shell aspects of the above type turn out to be important. The Einstein 
frame differs from the Jordan one by a field dependent conformal rescaling of the metric 
and a very nonlinear change of the Higgs field - radial variable of the Higgs multiplet. 
Quantization in these new variables has two major effects which are absent in the original 
Jordan frame. One effect is the change of the argument of the logarithms in quantum 
corrections, caused by conformal rescaling of mass parameters |^, This renders 

these logarithms nearly constant and not contributing to quantum corrections. Another 
effect is the change of coefficients of these logs [^]. One can show that not only Higgs but 
also Goldstone contributions get suppressed in the pre-logarithm coefficients. Moreover, the 
Einstein frame used in for the description of the chiral phase of the model automatically 
leads to the spherical parametrization of the Higgs multiplet, and this leads to an additional 
effect. The effective potential of the 0{N) multiplet, calculated in spherical coordinates, 
does not contain Goldstone contributions even off shell^^. Altogether, this leads to a big 
difference between the results in the Jordan and the Einstein frame - there are no large 
quantum corrections at high values of Mh when imposing the upper CMB bound on the 
Higgs range - a situation less sensitive to the CMB data, as claimed in 

This difference is interpreted in as an inalienable theoretical ambiguity of the quan- 
tization scheme. However, we are not ready to completely endorse this statement and 
prefer a pragmatic approach which is to try to reduce this ambiguity as much as possible. 
Since the cosmological parameters C, n-g, etc. are the functionals of the effective potential, 
this attempt can be based on the choice of a preferred parametrization which automatically 
renders this potential gauge-independent. Such a parametrization indeed exists - this is 
spherical coordinates for the Higgs multiplet. These coordinates disentangle the only local 
gauge invariant observable, the radial variable (p = [ip^-ip"-)^/'^ = ($^<I>)^/^. The effective 
potential for this variable is gauge independent, so that this solves both parametrization^^ 
and gauge independence issues, but of course leaves open the problem of the conformal 
frame. 

The price we will have to pay for such a resetting of the problem is that the theory 
becomes non-renormalizable even in the absence of gravitational interaction. Within such 
a setting it will make sense to include a two-loop RG improvement^^, establish error bars 
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On shell, y = 0, the values of the effective potential of course coincide, because of a vanishing Goldstone 



mass (2.6). 

^^Of course, one might ask the question of effective potential dependence on the choice of parametrization 
for Goldstone angles, but we do not even want to tread into the discussion of this conundrum. 

^'^ Two-loop RG improvement and associated with it error bars in |^ most likely exceed available precision 
for a number of reasons. Not the least of them is that the two-loop potential was taken from the calculations 
of in Cartesian coordinates, while its RG improvement in the chiral phase of inflating SM was based 
on beta functions derived in the spherical parametrization of the Higgs multiplet. Also, the conventional 
"Cartesian coordinates" RG flow of SM at low energy scales was matched with the chiral RG flow in 
spherical parametrization at the inflation stage. The consistency of this procedure requires checking both 



associated with the current indeterminacy of the SM and CMB data and the indeterminacy 
in the post-inflationary cosmological scenario. All this goes beyond the goals of the present 
paper, and will be considered in future publications. Here we are forced to reconcile with 
the gauge dependence of the obtained results. Fortunately, this dependence which enters 



through the renormalization function in (5.8) is rather weak. The value of in the 



Feynman gauge differs from that of the Landau gauge (used throughout the paper) by 
about 6%, which ultimately effects the value of spectral index in the third digit after the 
decimal point. 

Finally, in addition to a good match of the spectrum of cosmological perturbations 
with the CMB data our model also describes the mechanism of generating the cosmological 
background itself upon which these perturbations exist. This mechanism is based on the no- 
boundary or tunneling cosmological wavefunction [16, 17, 18, |2^] possibly prescribing 



initial conditions for inflation. Within the RG improvement the one-loop distribution 
function (0) IH for the no-boundary/tunneling state of the Universe is upgraded to 
p{ip) = exp ( =F r{ip)), where r{ip) is the (on-shell) value of the Euclidean effective action 
of the quasi-de Sitter instanton, weighting in the quantum ensemble the member with the 
initial inflaton value ip. The calculation of this action gives r{ip) = —2AT:'^Mp/V{ip) (see 
||ll[|) and the tunneling distribution function takes the form 

For the CMB range of Higgs masses this distribution features sharp peaks at the scale 
0g ~ ~ -647r2M|,/^o^7(io)/^ see Q for details. They exist only for A/ (to) < - 
which is exactly our case - and, of course, correspond to the peaks of the potential in Fig. 
|6| at t = 34 35. They have a very small quantum width Aip/ipQ ~ V^/C and, therefore, 
can be interpreted as generating initial conditions for inflation at the quantum scale ipo 
0, 13]. In the fuh CMB range of masses Mh their scale to exceeds the iV = 60 CMB 
formation scale tin, as it should for sake of chronological succession of initial conditions 
for inflation and the formation of CMB spectra (in fact, this is a direct consequence of 
the positive slope of V during inflation). Due to the high value of the exponentiated ratio 
i^^/A in ( [7.2| ) the probability peak is very distinctive even for the lower CMB value of 
Mh = 135.6 GeV. Moreover, for large Mh close to the upper CMB bound 185 GeV this 
peak gets separated from the non-perturbative domain of large over-Planckian scales due 
to a fast drop of y ~ A/^^ to zero shown on Fig. ^. This, in turn, follows from the fact 
that ^(t) grows much faster than A(t) when they both start approaching their Landau 
pole. Thus we get a complete tunneling prescription scenario of the quantum origin of the 
inflationary Universe [50| ]. 

To summarize, the inflation scenario driven by the SM Higgs boson with a strong non- 
minimal coupling to curvature looks very promising. This model supports the hypothesis 



at the analytical and numerical levels. 

Naively this differs by sign from the one-loop scale of inflation (3.15), but in fact this is a nontrivial 
equation for ipo because the right-hand side of this equality is a function of to. 



that an appropriately extended Standard Model can be a consistent quantum field theory 
all the way up to the quantum-gravity scale and perhaps explain the fundamentals of all 
major phenomena in early and late cosmology |51, Ultimately, it will be the strongly 
anticipated discovery of the Higgs particle at LHC and a more precise determination of 
the primordial spectral index Ug by the Planck satellite that might decide the fate of this 
model. 
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A. Renormalization of potential and non-minimal coupling terms — role 
of Higgs and Goldstone modes 

Here we justify the s-factor suppression mechanism due to back reaction of non-minimally 
coupled gravitons. For simplicity we consider only the Higgs multiplet-graviton sector of 
the model and the renormalization of the potential and non-minimal coupling terms of the 
classical action 

S = I d'xg''^ (u{m(9) - ^(V„<1.°)^ - , (A.1) 

U^'Jl±^, (A.2) 

For generality we work with the 0{N) scalar multiplet which for iV = 4 corresponds to 
the real representation of the SU(2) complex Higgs doublet. 

For renormalization of non-derivative terms of the action one expands it in quantum 
perturbations of the metric and scalar fields = [h^y, a) on the curved background with 
a constant background scalar field (^", = 99" + cr", 'Vf_if"' = 0. To cancel nonminimal 
derivatives in the kinetic term of the action (~ V^jV^y/i^'^) we pick up the background 
covariant DeWitt gauge-breaking term 

Sgb = -\ j d^xg'/^U{ip)(y>'h^, - lv,/i- ^(^aV.a'^)'. (A.3) 



The quadratic part of the total action then takes the form 



S2 + Sgb = -\ I d'xg'/'a^FABa' 



(A.4) 



with the following operator acting in the space of gravitons /i^jy and quantum Higgs mul- 
tiplet cr", the index A running over the range A = [fii^, a), 



AB 



a/3 



Sab + 



U 



UaUb ]n + Vab 



(A.5) 



where 



pa0 

111/ 



(A.6) 

(A.7) 
(A.8) 



Vab = m%{5ab - naUh) + mjjnarib - S,SabR- 

Here, the term —^-^naribO in the kinetic part of cr" originates from squaring the last 
term of the gauge, and ttiq = V l<~p and = V" are the Goldstone and Higgs masses 
given by Eqs. (^) and ( |2.7D above. They enter the potential term of the operator as 
coefficients of the projectors, respectively, transverse and longitudinal to the direction of 
the 0(A^)-multiplet vector of the background scalar field rf" = = ^J^°-Lp°' . 

The matrix coefficient of □ in F ab reads 



Cab 



(A.9) 



and its inverse equals 
^-lAB 



u 3u y/j-i^yaiB 



'-^riaUb 



-{Sab - nam) - SUaUb 



(A.IO) 



SO that the operator F = C takes the form 



6^ 



:R, 



F = ni + p- 

where the potential term consists of the two pieces 

P = P^+ pR. 

The first term is independent of the curvature and all its elements are 0(1/^) uniformly 
for the whole range of (p (or for all scales with 1 > s > 1/6.^) except the element 

1' 



(A.ll) 
(A.12) 



(P^)^ = -m^(5^ - n^Ub) - mjjs n^Ub + n^Ub O 



(A.13) 



It is dominated by the contribution of the Goldstone mass, and we exphcitly present the 
term showing the suppression of the Higgs mass contribution by the factor s. 
The second term is hnear in the curvature and reads 



(A.14) 



where 



'R)b 



Now we calculate the one-loop effective action 



R. 



(A.15) 



^i-loop = 1 Tr In [ ni + P - -^i? 



(A.16) 



in the approximation linear in spacetime curvature. Its divergent part is given by the 
trace of the second Schwinger-DeWitt coefficient tr 02 [^. In this approximation it is just 
one half of the square of the full matrix potential term ( A.12| )-^'^. Since the -P^ part of this 
term commutes with the rest of the operator (A. 11) in the 1/^-approximation (in view of 
the matrix structure and in view of the assumed constancy of the background V^tp = 0), 
it can be considered as a mass matrix 



(A.17) 



so that 



(A.18) 



2 2 

and the one-loop action in dimensional regularization, (i ^ 4, takes the form explicitly 
featuring the Coleman-Weinberg structure of the effective potential, 



1 



327r2 I 2 



-|- In 47r tr 0,2 



1 /,v2^2/, 3 
tr (M ) I In 



647r2 
1 

' 32^ 



2 2 



tr PijM^( In 



M 



1 + 



(A.19) 
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Because of the constant background scalar field, the Faddeev-Popov ghost contribution is (^-independent 



and does not contribute to the renormalization of the sector (A.l) of the model. 

^"With a non-constant background V^^p 7^ the operator (A_.ll) acquires terms linear in derivatives. 



which can be absorbed into the generalized d'Alembertian by the redifinition of the covariant derivative, 
= -I- 0(\7f^(fi). Then, extra terms in 0,2 are given by the square of the curvature associated 
with the commutator of [1)^,73^] and D^P. They are responsible for the renormalization of the kinetic 
term and curvature squared terms in the effective action, which we disregard here. 



Here, indeed, the second line represents the Coleman- Weinberg potential and the third line 
gives the logarithmic renormalization of the non-minimal coupling, because after calculat- 
ing the trace it takes the form U^~^°°^{ip)R. Since M ~ 99^, the logarithms of the mass 
matrix yield, up to subleading corrections, a unit matrix llnc^^ and the traces equal 



647r2 ^ ' 327r2 



m^{N - 1) + s'^m% 



+ 



^ aV 

~ 1287r2 

1 ^ ,',2 1 



;6 + 18s2)+oQ^ (A.20) 



trPfiM = —^R 



327r2 " 327r2 



mliN - 1) + sV^J C + 0(C°) 
3A 

= v'R^{i+s')c+o{e). (A.21) 



These equations show the suppression of the Higgs-mass contribution by the factors of s 
and the origin of 6 A term of the anomalous scaling in Eq.( |3.2D (cf. Eq.( |2.8| )). This term is 
entirely due to the contribution of Goldstone modes. Similarly, they explain the origin of 
the dominant part of the constant C ~ ^ in ( |3.3[ ) , given by (^) . 

These equations also underlie the beta functions ( |4.8D and ( [4.9D in which we retained 
only the leading terms in ^ ^ 1, nontrivially depending on s-factors. One can check that 
subleading terms of ( [A.20 ) and ( [A.21 ) above are uniformly damped by inverse powers of 



^ in the full range scales 1 > s > 1/6^. They originate mainly due to graviton loops 
suppressed by powers of effective Planck mass l/(Mp -|- ^v?^). The remaining part of the 
beta functions ([l.S|)- ([4.9|) is generated from the heavy vector boson and quark sectors by 
adding to the mass matrix their contribution diag [m^^, m^, ] and by including 
the anomalous dimension terms caused by field renormalization. 
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